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Abstract 

High-temperature superconductivity in both the copper-oxide and the iron-pnictide/chalcogenide systems 
occurs in close proximity to antiferromagnetically ordered states. Neutron scattering has been an essential 
technique for characterizing the spin correlations in the antiferromagnetic phases and for demonstrating 
how the spin fluctuations persist in the superconductors. While the nature of the spin correlations in the 
superconductors remains controversial, the neutron scattering measurements of magnetic excitations over 
broad ranges of energy and momentum transfer provide important constraints on the theoretical options. 
We present an overview of the neutron scattering work on high-temperature superconductors and discuss 
some of the outstanding issues. 



1. Introduction 

Looking at temperature vs. composition phase 
diagrams of high-temperature superconductors, in- 
cluding both cuprate and iron-pnictide/chalcognide 
systems, one finds antiferromagnetic (AF) order in 
close proximity to superconductivity [THl]. This 
close association, together with experimental evi- 
dence for AF correlations in superconducting sam- 
ples, has led many theorists to a common belief 
that AF fluctuations play an important role in the 
electron-pairing mechanism that underlies super- 
conductivity [3H5] • The differences among theoreti- 
cal perspectives only begin to appear when one con- 
siders the speciflcs of how the AF correlations im- 
pact pairing. To understand the source of these dis- 
agreements, one must step back and recognize that 
there is no consensus on how to describe the inter- 
action between charge carriers and spin excitations 
in a metallic conducting system with strong AF cor- 
relations, especially when the system is close to a 
transition to a correlated (Mott) insulator state. In 
the absence of a common view on how to frame the 
problem, it should not be surprising that there is a 
lack of consensus on the approach to a solution. 

Neutron scattering is the preeminent technique 
for studying AF spin fluctuations in solids. One 
can, of course, also obtain important information 
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on local hyperfine fields, susceptibilities, and slow 
fiuctuations from techniques such as nuclear mag- 
netic resonance (NMR) and muon spin rotation, 
and magnetic order can be probed with resonant 
x-ray scattering. The defining feature of AF fluctu- 
ations in high-temperature superconductors, how- 
ever, is their remarkably high energy scale. Mag- 
netic excitations in these systems extend up to en- 
ergies of several hundreds of meV, which easily ex- 
ceeds the maximum energy of phonon excitations 
involved in the traditional mechanisms of supercon- 
ductivity. This observation makes highly energetic 
AF fluctuations a primary suspect for mediating 
the unconventional high-temperature superconduc- 
tivity. 

To characterize the AF correlations in these and 
other strongly-correlated metallic systems, it is 
therefore crucial to cover the energy range from 
fractions to hundreds of meV and to probe all 
of reciprocal space. Developments in spectrome- 
ters and sources over the last two decades have 
greatly improved the efficiency of such measure- 
ments. The throughput of triple-axis spectrometers 
has been enhanced through the use of multiple ana- 
lyzers and detectors, beginning with developments 
such as RITA [B| (first at Ris0 and continuing at 
SINQ) and SPINS [7] and progressing to MACS [S] 
and BT7 at the NIST Center for Neutron Re- 
search; of course, enhanced instruments are avail- 
able at many other facilities, including the Insti- 

January 25, 2013 



tut Laue Langevin, FRM-II, and the Laboratoire 
Leon Brillouin. While most triple-axis instruments 
are optimized for excitations in the cold to thermal 
neutron range {E < 30 meV) time-of-flight spec- 
trometers at spallation sources have enabled mea- 
surements with cpithermal neutrons probing exci- 
tations up to ^ 1 eV, and covering a very large 
phase space. This began with the seminal MAPS 
spectrometer [TOj (and later MERLIN [IT]) at ISIS, 
and has become common place with instruments 
such as ARCS ^ and SEQUOIA [TT at the Spal- 
lation Neutron Source (SNS) and 4SEAS0NS [H] 
at the Japanese Proton Accelerator Research Com- 
plex. Very soon, the SNS spectrometer HYSPEC 
will have the capability of distinguishing magnetic 
from nuclear scattering through neutron-spin polar- 
ization analysis [151 , a technique that has long been 
used to advantage on triple-axis instruments |16j . 

Ideally, experimentalists should be able to inter- 
pret their results by comparing with existing theo- 
retical expressions or to simply measure the mag- 
netic response in a superconductor family as a func- 
tion of doping and temperature, and then let the 
theorists interpret the results. In practice, how- 
ever, things are not so simple. The neutron scat- 
tering cross section is relatively weak, so that large 
crystals (> 1 cm'^) or assemblies of crystals are 
typically required for each composition to be stud- 
ied. In some cases, it has taken decades for crystal 
growth technology and know-how to evolve to the 
point that suitable crystals are available. Another 
aspect involves interpretation of the measurements. 
Frequently, there is a disconnect between theoreti- 
cal interpretations and experimental results arising 
from various assumptions and idealizations adopted 
both in analyzing the data and applying the theo- 
retical model. To present results in a useful fashion, 
it is generally necessary to parametrize the data by 
fitting with a model. The models used, or the words 
used to describe them, are often based on specific 
theoretical perspectives. Without a general and ac- 
cepted theoretical description of itinerant antifer- 
romagnetism, there is no universally accepted and 
unbiased language for describing measurements of 
spin fluctuations in high-temperature superconduc- 
tors. 

In the following sections, we will briefly summa- 
rize what we consider to be the important results 
from neutron scattering studies on antiferromag- 
netism in hole-doped cuprates and in iron-based 
superconductors. We will also point to some of the 
open questions. Obviously, our choices and presen- 



tation reflect our own biases. The coverage of ex- 
periments and theory is necessarily incomplete, and 
we refer the interested reader to more extensive re- 
view articles on cuprates [H [T7H?7] and iron-based 
systems [IUHHSI]. 

2. Notes on neutron scattering 

Neutron scattering measures the product of the 
dynamical spin structure factor 5(Q,w), which is 
the Fourier transform of the spin-spin correlation 
function, and the square of the magnetic form fac- 
tor F(Q), which is the Fourier transform of the den- 
sity of the electronic magnetization cloud associated 
with each spin, normalized to one at Q = 0. Here 
/iQ is the neutron momentum transfer and —fku is 
the energy transfer. The dynamical structure fac- 
tor describes the cooperative behavior of electronic 
spin variables, including spin order and excitations, 
whereas the form factor relates it to the behavior of 
the magnetization density in the crystal. This latter 
is the quantity which interacts with neutron's mag- 
netic moment and is probed in experiment. The 
magnetic form factor is determined by the Wannier 
functions of magnetic electrons, which could be ob- 
tained from first-principles calculations. A predic- 
tion for the dynamical structure factor, 5(Q,(u;), 
can in principle be obtained from the theoretical 
analysis of the model spin Hamiltonian that de- 
scribes the system. For an ordered antiferromag- 
net, the excitations are spin waves that can be cal- 
culated in perturbation theory, where one assumes 
that the fluctuations are small compared to the or- 
dered moment. As we will discuss, the fluctuations 
become more important when one reduces the size 
of the spin and the dimensionality. In the case of 
the weakly-coupled planes with spin S ^ that 
occur in cuprates, spin-wave theory is not a well- 
controlled approximation, though it remains a use- 
ful description. 

5(Q,a;) obeys a number of useful sum rules, 
which are helpful in analyzing the experimental 
neutron scattering data [7]. In particular, for a sys- 
tem of identical spins S, one has 

cf dio f dQ5(Q,a;) = 5(5 + 1), (1) 

where C = uo/(27r)'^ and vq is the unit cell volume. 
It is possible to convert S{Q,uj) to the imaginary 
part of a generalized spin susceptibility with the 
formula 

x"(Q,^) = 7r(l-e-''"/^-^)5(Q,c.). (2) 



By integrating x"{Qt^) over Q, one obtains the 
local susceptibility x"{'^)- 

3. Cuprates 

3.1. Parent insulators 

The parent materials of the cuprate supercon- 
ductors are correlated insulators. Taking La2Cu04 
as an example, it has a charge excitation gap of 
~ 2 eV [TB], so that the only low energy excita- 
tions involve spin fluctuations. These can be de- 
scribed by an effective spin Hamiltonian, where the 
nearest-neighbor Cu spins within the planes are 
coupled by a superexchange energy J that is quite 
large, while the effective coupling between planes 
is extremely weak, so that there are strong two- 
dimensional (2D) AF spin correlations at tempera- 
tures far above magnetic ordering temperature 
[33]. 

The Cu moments within the Cu02 planes or- 
der antiferroniagnetically below a Neel temperature 
T/v — 325 K in stoichiometric La2Cu04 [31] ■ Neu- 
tron diffraction studies have shown that the ordered 
moment is quite small, ~ 0.6^b |35j . This cor- 
responds to an ordered spin value of (S) w 0.27 
(assuming g « 2.2 for the spectroscopic Lande 
t;— factor of Cu [33]), which is even smaller than the 
value of 0.30 predicted by spin- wave theory upon 
accounting for zero-point motion |37H39j . Similar 
results have been obtained for other cuprate anti- 
ferromagnets [22J- Consequently, the contribution 
of the static spin order, which is measured by elas- 
tic Bragg peaks (delta- functions in energy), (S")^, 
accounts for only ~ 10% of the total spin spectral 
weight, S{S+1) = 3/4, in the sum rule for 5" = 1/2. 
Hence, more than 90% of the spectrum is inelastic. 

Figure [T] shows a recent measurement of the spin 
excitations in antiferromagnctic La2Cu04 by Head- 
ings et al. [1^. The authors find that both the 
dispersion and the wave-vector dependence of the 
intensity are described rather well by spin- wave the- 
ory expressions. The absolute intensity, however, 
appears much lower than predicted by the linear 
spin-wave theory, requiring a downward renormal- 
ization by a factor of Zd = 0.4 ± 0.04. This is 
somewhat smaller than the value Zd ~ 0.6 that is 
predicted from quantum corrections to linear spin 
waves [IT] . 

A fit to the dispersion in La2Cu04 yields J = 
143 ± 2 meV, but also requires longer-range ex- 
change couplings — second and third neighbor cou- 
plings J' and J", which are relatively weak, and a 
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Figure 1: (a) Spin wave dispersion in La2Cu04 measured 
by inelastic neutron scattering, plotted for high-symmetry 
directions, as indicated in the inset of (c). (c) Measured 
single magnon intensity vs. wave vector. In both panels, the 
lines represent fits with a spin-wave model. Reprinted figure 
with permission from Headings et al. |40| . @ 2010 American 
Physical Society. 



significant 4-spin cyclic exchange term Jj, « 0.4J. 
The overall bandwidth of the magnetic spectrum is 
~ 2 J. At the highest energies, there are some mod- 
est deviations in the spectral shape relative to the 
single-mode spectrum predicted by spin-wave the- 
ory. These seem to be consistent with corrections 
obtained in quantum Monte Carlo calculations [32] ■ 
Similar findings were recently reported for several 
other cuprates by Dalla Piazza et al. [H] , where the 
spin excitation spectra can be well described by a 
perturbative (up to a second order) treatment of an 
effective 1-band Hubbard model with nearest- and 
next-nearest neighbor hopping, or by a spin-wave 
treatment of the local spin Heisenberg Hamiltonian 
with extended-range interactions J, J', J" and Jc- 
Quantitative estimates for J require that one take 
account of the bridging O atoms when evaluating 
the Cu to Cu hopping [33]. The 4-spin cychc ex- 
change can be obtained at the same order of ap- 
proximation in such a multi-orbital model [IS]. 

Evaluating the sum rule for La2Cu04 without 
any model assumptions for iS(Q,a;), the result is 
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Figure 2: (a) Constant magnitude surface of the magnetic 
form i^{Q) calculated from the Wannier function for a Cu 



3d^ 



orbital hybridized with its neighbors in Sr2Cu03. 



(b) Similar plot for an ionic Cu 3d^2_j,2 state, (c), (d) Com- 
parison of cuts of the magnetic form factors shown in (a) 
(solid blue line) and (b) (dashed red line) along two symme- 
try directions. From Walters et al. ^48 . 
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Figure 3; Magnetic dispersion relation along 

Q = (0.5 -I- h, 0.5, 0) in various cuprates, correspond- 
ing to wave vectors parallel to the Cu-O bonds. Data 
are for Lai.875Bao.i25Cu04,[50] Lai saSrn iBCuOd.l51| 
YBa2Cu306.6,152j YBa2Cu306.5,[521 iHl and 
Bi9Sr2CaCu9 0K+.s.l55| The energy is scaled by J for 
the AF parent material. [22ll56] From Fujita et al. [T7] . 



~ 60% of the expected result for S" = i gO]. Al- 
though part of the missing weight could be in multi- 
magnon excitations that are out of the measure- 
ment range, a significant missing weight may result 
from not taking proper account of the distribution 
of spin density. Shamoto et al. 46j] have shown 
that accounting for the Zd^2_y2 orbital anisotropy 
of the Cu^+ magnetic form factor is absolutely es- 
sential to understand the magnetic Bragg diffrac- 
tion in YBa2Cu306.i5; it is also important for the 
analysis of the spin dynamics [40l |47]. Although 
using the anisotropic ionic magnetic form factor of 
Cu^"*" is much better than using a spherical form 
factor of the dipole approximation, it still neglects 
the effect of covalency, i.e. charge transfer to the 
oxygen neigbors, which turns out to be very signif- 
icant in the cuprates. 

Walters et al. [48] did a careful study of the form 
factor in the quasi-lD antiferromagnet Sr2Cu03. 
Local structure of the planar Cu-O square plaque- 
ttes in this material is essentially identical to that 
in La2Cu04. Making use of a precise theoretical re- 
sult for the excitation spectrum, they demonstrated 
that a good fit to the data requires a form fac- 
tor that takes account of hybridization between the 
half-filled Cu 3dx2_y2 orbital and the ligand O 2po- 
orbitals, as given by a density functional calcula- 
tion. The hybridization causes the spin density to 
be extended in real space, resulting in a more rapid 
fall off in reciprocal space compared to a simple 



Cu^+ form factor, as illustrated in Fig. [2j Smaller 
values of magnetic form factor at the wave vectors 
where the measurement is performed lead to the 
suppression of magnetic intensity, which could be 
as large as a factor of two or more [48]. Finally, we 
note that a study of covalent NMR shifts in by Wal- 
stedt and Cheong 49J found that barely 2/3 of the 
spin density in La2Cu04 resides on the copper sites, 
in excellent agreement with the Sr2Cu03 neutron 
data of Walters et al. [48j . 

3.2. Hole-doped superconductors 

The antiferromagnetic insulator state is relatively 
well understood. To get superconductivity, it is 
necessary to dope holes into the Cu02 planes. 
What happens to the spin fluctuations? 

Figure [3] summarizes the effective magnetic dis- 
persion found in underdoped to optimally-doped 
cuprates. The energy scale is normalized to J 
of the undoped parent materials. (Note that J 
can vary by 30% among different cuprate families 
|56L I57j .) There are two important observations to 
make regarding this hour-glass-like spectrum: 1) 
the scale of the high-energy excitations seems to be 
J, characteristic of the correlated insulator state, 
and 2) the low-energy excitations are incommensu- 
rate. Concerning point (1), the energy of the neck 
of the hour glass, i?cross, decreases towards zero as 
one reduces the doping and approaches the anti- 
ferromagnetic state [T7]. Also, as we will discuss 
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shortly, the spectral weight of the magnetic excita- 
tions for E > E'cross is observed to evolve smoothly 
from the AF state. 

For E < -Ecross, samples near optimal doping and 
above develop a gap in the spin excitations and a 
resonance peak when the temperature drops below 
the superconducting transition, Tc- An example of 
the change in the magnetic response at the AF wave 
vector for nearly-optimally-doped YBa2Cu306+a; is 
shown in Fig. g] Here the magnetic signal is pre- 
sented in terms of the imaginary part of the dy- 
namical spin susceptibility, x"(Q,'^)- The energy 
of the resonance, Er, varies with doping and among 
different cuprate families. From measurements on 
YBa2Cu306-i-a; and Bi2Sr2CaCu208+5, Sidis et al. 
[59] found that Er « 5.3fcTc. Alternatively, Yu et 
al. [BU] compared with the energy scale A of the 
d-wave superconducting gap determined in photoe- 
mission studies and found Er/2A w 0.64 for a large 
variety of materials. 

How does one explain the temperature- 
dependent spin gap and resonance peak? One 
way is to assume that the same electrons that go 
superconducting are also responsible for the mag- 
netic response. In this case, the magnetic response 
is presumably due to scattering the conduction 
electrons across the Fermi surface. If the Fermi 
surface is "nested", meaning that it has parallel 
portions separated approximately by the AF wave 
vector, then one can obtain a substantial spin 
response [ST], as occurs in the spin-density- wave 
state of metallic chromium [H5]. The opening of 
the superconducting gap removes electronic states 
that can contribute to x", thus resulting in the 
spin gap. The resonance peak is a consequence of 
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Figure 4: Measurements of x"{QaFi'^) below and 
above Tc demonstrating the resonance and spin gap in 
YBa2Cu306.92; from Bourges et al. |58l . 



BCS coherence factors plus a superconducting gap 
that changes sign between the nested portions of 
the Fermi surface [3j. This is consistent with the 
d-wave gap of the cuprates [53] . 

There have been many theoretical calculations of 
the resonance, as well as the downwardly dispers- 
ing excitations below it They can describe 
rather well the resonance in YBa2Cu306+£c, which 
is commensurate with the AF wave vector, and they 
can provide some qualitative agreement with the 
downwardly-dispersing excitations; however, there 
are also significant discrepancies. For T <^Tc and 
huj ^ 2A, the magnetic response should be deter- 
mined by the positions of the nodes of the super- 
conducting gap ^64jj but experimentally the low- 
energy incommensurate scattering is rotated 45° 
from the predicted direction pS], and there is no 
sign of the predicted signal even at the lowest ener- 
gies |66j . Furthermore, the resonance in optimally 
and over-doped La2_a;Sr2;Cu04 occurs at incom- 
mensurate wave vectors [501 [HZ], rather than the 
commensurate position typically predicted by the 
calculations [55]. 

An alternative explanation for the low-energy in- 
commensurate spin excitations attributes them to 
an electronically heterogeneous state involving al- 
ternating charge and spin stripes [531 [BS] [53] . Static 
ordering of charge and spin stripes is observed in 
La2_a:Ba2,Cu04 [7^ and related systems [l7| in 
which a subtle structural transition breaks the ro- 
tational symmetry of the planar Cu-0 bonds, al- 
lowing a unique orientation of the stripes within 
each layer. At low energies, spin fluctuations rise 
out of the incommensurate magnetic superlattice 
peaks, and at higher energies they exhibit the hour- 
glass spectrum of other cuprates, as shown for 
La2-a:Ba2,Cu04 with x = 1/8 in Fig. [sj Although 
stripe order tends to compete with 3D supercon- 
ducting phase order, it coexists with 2D supercon- 
ducting correlations [TTJlTj] and even layered phase- 
decoupled superconductivity [73]. It has been pro- 
posed that the coexisting spin and superconduct- 
ing orders are intertwined in a pair-density-wave 
(PDW) state [71]. In such a state, the Cu moments 
providing the spin response are distinct from the 
doped holes that form the PDW superfluid. 

It is worth noting that at higher temperatures, 
where the stripe order is lost, the spin excitations 
remain incommensurate [751 176] and are quite simi- 
lar to those in optimally-doped La2_a:Sr2:Cu04 [77]. 
The anisotropic bulk susceptibility in paramagnetic 
La2_2;Ba2,Cu04 is consistent with local moment be- 



5 



havior [36 . In YBa2Cu306+x, the modulation of 
the magnetic response as a function of momentum 
transfer perpendicular to the Cu02 bilayers indi- 
cates that the relevant spins are on Cu sites [75] 
and not on O, where the doped holes are located 

Another place to look for the impact of doped 
holes is at high excitation energies. A num- 
ber of groups have used time-of-flight measure- 
ments to characterize the spin fluctuations over 
a broad energy range in YBa2Cu306+x [HOI I5T] . 
La2-xSra;Cu04 EH EMSl, La2-a;Ba:rCu04 [SHI 
[76] . and Bi2Sr2CaCu208+5 [55]. Analyzing the lo- 
cal susceptibility, x"i^)^ Stock et al. [54] noticed 
that the strength of the magnetic response falls off 
dramatically and systematically above a doping- 
dependent energy scale. To quantify this behav- 
ior, they evaluated the energy at which x"(w) falls 
below half of that for an undoped antifcrromag- 
net; their results are shown in Fig. [5] Remark- 
ably, the identified energy scale corresponds very 
well with the pseudogap energy determined from 
electronic spectroscopies [SS]. The resulting pic- 
ture is that spin fluctuations retain the strength of 
a correlated insulator below the pseudogap energy, 
where charge excitations are reduced, while mag- 
netic weight is strongly depressed above the pseu- 
dogap, where charge excitations are stronger and 
better defined. 

A recent resonant inelastic x-ray scattering 
(RIXS) study on YBa2Cu306+2; and closely related 
compounds fSS' appears to confiict with these re- 
sults. First note that because of limitations on mo- 
mentum transfer and energy resolution, the RIXS 
measurements are limited to Q = (h, 0, 0) with h < 
0.45 reciprocal lattice units and energies greater 
than 130 mcV. The excitations observed in an in- 
sulating sample are consistent with expectations 
for antiferromagnetic spin waves and appear con- 
sistent with neutron results. This is reasonable, as 
the presence of antiferromagnetic order ensures that 
the dispersion must be the same about both ferro- 
magnetic and antiferromagnetic zone centers. Mea- 
surements on superconducting samples up to opti- 
mal doping, however, find excitations that soften no 
more than 10% with doping, with negligible change 
in integrated intensity. Such results are quite differ- 
ent from the neutron scattering results, especially 
those shown in Fig. [5j they also conflict with 2- 
magnon Raman scattering results from Sugai et al. 
[56l . Of course, in the absence of antiferromagnetic 
order there is no constraint that the excitations, 
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Figure 5: Large symbols: estimates of doping-dependent 
energy scale at which magnetic spectral weight falls below 
half that of the antiferrromagnetic state based on inelastic 
neutron scattering studies of various cuprates. Small gray 
symbols: pseudogap energy from various electronic spectro- 
scopies as summarized in 1851 . Reprinted figure with per- 
mission from Stock et al. '54', © 2008 American Physical 
Society. 



especially including the effects of damping, be the 
same in both the ferromagnetic and antiferromag- 
netic zones. Now, the RIXS cross section includes 
both charge and spin excitations, and it is not yet 
firmly established how the charge and spin channels 
will interact in measurements of metallic samples. 
Even assuming that a proper identification of mag- 
netic excitations has been made, there is no phys- 
ical justification for the assumption made in [86] 
that the dispersion measured by RIXS near Q = 
is the same as that near the antiferromagnetic wave 
vector, where the neutron measurements have been 
done. The neutron cross section is very well un- 
derstood, and the indications of strong damping of 
large Q magnetic excitations are clear. As the large 
Q excitations are the ones relevant to theory, the 
significance of the RIXS measurements is not clear. 

Implicit in Fig. [5] is the fact that AF excitations 
die away with over doping. An explicit demonstra- 
tion of this is given in Fig. |6] from |83j ; related re- 
sults for La2_a;Sra;Cu04 with x — 0.22 have been 
reported by Lipscombe et al. [82 . Thus, as one 
dopes far away from the AF insulator state, the AF 
spin fluctuations disappear. This occurs despite the 
fact that the Fermi surface becomes better nested 
with overdoping in La2-i:Sr3,Cu04 87J. 
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Figure 6: Comparison of magnetic intensity vs. energy in 
over-doped La2-a;Sra;Cu04 with x = 0.25 (red diamonds) 
and X = 0.30 (black circles) with La2-a;Baa;Cu04 with x = 
0.125 (blue squares). Fom Wakimoto et al. 83 , © 2007 
American Physical Society. 



To summarize, experiments indicate that spin 
fluctuations in the cuprates evolve continuously 
from the AF insulator state. Hole doping causes 
changes at low energy — incommensurate disper- 
sion, spin gap, and resonance — and at high energy, 
where the spectral weight drops above the pseu- 
dogap energy; the AF fluctuations disappear with 
overdoping, where the superconductivity also goes 
away. Model calculations in which the doped car- 
riers are assumed to provide the magnetic response 
provide an appealingly direct connection to the su- 
perconductivity; however, they have a challenge 
to explain the doping dependence of the magnetic 
weight and the role of the superexchange energy 
J. Alternatively, phase separation of the holes into 
charge stripes provides a natural way for locally- 
AF spin correlations among Cu moments to survive; 
the challenge here is to explain, from first principles, 
the transport properties and the superconductivity. 
Empirically, the differences in the electronic proper- 
ties of cuprates with stripe order and those without 
is extremely subtle Models of pairing [ST] 

and superconductivity [22] in striped systems have 
been proposed, but those models tend to involve ad 
hoc features. 

Before leaving the cuprates, we should touch on 
the topic of "intra-unit-cell" magnetic order. As 
reviewed by Bourges and Sidis [93 , a proposal of 
loop-current order by Varma 94 has motivated po- 
larized neutron diffraction searches for a potential 
hidden order. Varma's current loops are between 
a Cu atom and nearest-neighbor oxygens, so that 



the scattering should appear at Bragg peak posi- 
tions; however, the size of the current loops implies 
a form factor that falls off rapidly with Q. In fact, 
the spatial extent of these current loops is similar to 
that of the strongly extended Wannier functions of 
magnetic electrons in cuprates, revealed by the co- 
valent magnetic form factor [JS]. Hence, the extent 
of the form factor of the loop contribution should be 
similar to that of the spin magnetic moment con- 
tribution, but the direction of the loop magnetic 
moment must be determined by its orbital origin. 

Neutron studies have found evidence for a change 
in polarized neutron scattering for certain low-Q 
peaks in several cuprate families below tempera- 
tures consistent with the onset of the pseudogap. 
In addition, Li et al. |95j have reported a high- 
energy, dispersionless, apparently-magnetic mode 
in HgBa2Cu06, that has been interpreted as an 
Ising-like excitation of Varma's state [96] . We note 
that Lederer and Kivelson [97^ have done an analy- 
sis of the loop current model and shown that, in 
combination with the neutron results, it implies 
finite magnetic fields at various lattice sites that 
should be detectable by nuclear magnetic resonance 
(NMR). Definitive NMR tests have not yet been re- 
ported. 

We are uncertain what to think of the measure- 
ments of intra-unit-cell magnetic order. The mea- 
surements themselves are challenging to perform, 
as they require identifying a small magnetic signal 
on top of a substantial nuclear Bragg intensity. The 
observed polarization [98^ of this additional contri- 
bution is at odds with the original model of current 
loops in the a6-plane :94 . The situation would be 
easier to judge if there were other systems where 
this sort of magnetism has been observed. (A re- 
cent experimental claim [99] of orbital currents in 
antiferromagnetic CuO has been explained away in 
terms birefringence effects of the monoclinic lattice 
[100) .) In contrast, stripe order has been observed 
in a number of transition-metal oxides [101| . and, 
of course, there are examples of metallic antiferro- 
magnets, such as Cr, where Fermi-surface nesting is 
important [62'. If orbital currents are real, how do 
they interact with the spin response that has been 
clearly identified? 

4. Fe-based superconductors 

Nature of the magnetism 
The crystal structure of iron pnictide and chalco- 
genide superconductors features a square-lattice ar- 
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rangement of magnetic ions, similar to that in 
cuprates [2l l28lfT02Ul05) . It consists of a continuous 
stacking of square-lattice layers of iron atoms, each 
sandwiched between the two half-density layers of 
bonding pnictogen or chalcogen anions. These an- 
ions, which tetrahedrally coordinate the Fe sites, 
occupy alternate checkerboard positions above and 
below the Fe layer, so that the resulting unit cell 
contains two formula units. The sparsely spaced 
layers are only weakly held together in this quasi- 
two-dimensional structure, which facilitates the re- 
duced structural symmetry. 

Unlike in the cuprates, where bridging oxygen an- 
ions connect the nearest-neighbor copper sites thus 
mediating the dominant nearest-neighbor superex- 
change interaction, bonding anions in the structure 
of iron-based superconductors connect the sites on 
the diagonal of the square plaquette. This suggests 
that the next-nearest neighbor, diagonal coupling 
might be the dominant superexchange interaction 
in the Heisenberg Hamiltonian of the local-spin pic- 
ture for iron pnictides and chalcogenides |1061 1107| . 
The resulting J1-J2 Heisenberg model exhibits frus- 
tration and has a complex phase diagram, includ- 
ing disordered spin-liquid and different Neel states 
[108 ' . In addition, the nearest Fe-Fe distance is ex- 
traordinarily short, only « 2.8 A, suggesting that 
direct orbital overlaps could also be important, and 
so might be further-neighbor couplings. Extending 
the interactions to include coupling of second neigh- 
bors along the side of the square leads to a J1-J2-J3 
model on a square lattice jl09j . Understanding dif- 
ferent non-classical frustrated phases and the corre- 
sponding spin excitations in this model presents a 
major challenge. Thus, if we consider iron-based 
materials in terms of effective local-spin models, 
the situation is substantially more complex and less 
well understood than in the cuprates. 

The most important distinction between the 
cuprates and the iron pnictides and chalcogenides, 
however, stems from the 3dP electronic configura- 
tion of the Fe atom, compared to the 3(P config- 
uration of Cu |110| . Hence, while a single dx2-y2 
orbital is involved in describing the band structure 
of the cuprates, there are at least four active d- 
orbitals, in addition to the anion bonding p-orbitals, 
contributing to the unoccupied, or partially oc- 
cupied band structure in the iron systems. In a 
strongly-correlated, localized-electron ionic picture, 
four d-orbitals would result in an 5' = 2 spin system, 
if Hund's rule is obeyed. Experiment, however, re- 
veals maximum spins of = 3/2 in the most mag- 



netic system, FeTe 1112] , with smaller effective 
spins in ferropnictides pll28j. Half- integer spins im- 
mediately rule out the purely local-spin ionic pic- 
ture with localized d-electrons, as only integer-spin 
states are allowed in this case. This, of course, 
agrees with the metallic conductivity in iron pnic- 
tides and chalcogenides observed in experiment. 

It is perhaps worth recalling that the effective lo- 
cal spin is revealed in neutron scattering through 
the total scattering intensity. As given by Eq. ([T]), 
the integrated spectral weight corresponds to 5(5-1- 
1) = [/icff/(5/^B)]^, where g is the Lande factor. 
This defines the fiuctuating instantaneous effec- 
tive moment, /loff, whereas the ordered static mo- 
ment measured by Bragg diffraction in the mag- 
netically ordered state, (/i) — gfiB{S), is always 
smaller. Here (5) is the ground-state value of the 
spin operator S, and g^siS) < gfJ-sS < /ieft- In 
low-dimensional and/or frustrated systems, such as 
cuprates and iron pnictides and chalcogenides, (S) 
could be much smaller than the actual local spin S. 

In the band theory description, even the num- 
ber of electrons (holes) precludes iron-based ma- 
terials from being Mott insulators and renders 
them semiconductors or semimetals with electron 
and hole pockets. Understanding the properties 
of these systems requires a multi-band descrip- 
tion from the outset [113J. Although many treat- 
ments of the itinerant-electron limit consider only 
two or three t2g derived bands formed by dxz, dyz 
and dxy orbitals hybridized with the p-orbitals of 
the anion [TT^ 1115) . some researchers suggest that 
such treatment is insufficient and that more bands 
should be included |116j . Theoretical results includ- 
ing all five d-orbital derived bands have been used 
to describe the neutron scattering data |117H119] . 
The weakness of the itinerant-electron description 
is that it only accounts for the on-site interactions 
by either applying perturbation theory [ 120.^ , which 
fails when the interactions are substantial, or within 
the mean- field random phase approximation fll7]. 
which is not well controlled. Perhaps the most ade- 
quate way of treating the magnetic response of iron 
pnictides and chalcogenides theoretically, is by em- 
ploying the dynamical mean field theory (DMFT), 
or a combination of DMFT with the density func- 
tional theory (DFT) P^THHS] . DMFT is designed 
to account for the on-site correlation, which in pnic- 
tides and chalcogenides is controlled to a large ex- 
tent by the spin-dependent Hund's rule coupling 

m- 

To date, there have been numerous studies of 
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Figure 7: Magnetic excitations from BaFe2As2, modeled with a Heisenberg spin Hamiltonian by varying exchange parameters 
Jla, Jib} S'lid J2, and applying an anisotropic damping. Reprinted figure with permission from Harriger et al. [126| . © 2011 
American Physical Society. 



the magnetic excitations in the antiferromagnetic 
members of various Fe-based superconductor fam- 
ihes H [Ml EH [USHIin]. In many papers, there 
is a focus on determining whether the excita- 
tions are better described in terms of the spin- 
wave spectrum derived from a spin-only Hamilto- 
nian or of a weakly-coupling itinerant-calculation 
involving Fermi-surface nesting. We argue that 
this is an ill-posed issue. Calculations of a 
typical antifcrromagnetically-ordered state within 
the local-spin density approximation indicate that 
there is spin polarization of states at high binding 
energies (~ \eV) |130| . This method does not work 
in the paramagnetic state; however, attempts to 
calculate the phonon energies only find agreement 
with experiment when lattice constants are used 
that are consistent with those obtained in calcu- 
lations with magnetic order |131j . (Spin-polarized 
ions are bigger than unpolarized ions.) The elec- 
tronic interactions, dominated by Hund's rule cou- 
pling, tend to be intermediate in strength |1141I132] . 
so neither weak-coupling nor localized-moment ap- 
proaches are accurate. 

Nevertheless, it is generally reasonable to use 
a spin-wave model based on a suitable spin-only 
Hamiltonian to parametrize measurements. Fig- 
ure[7]shows the spin excitation spectra measured by 
Harriger et al. |126| in the AF phase of BaFe2As2. 
As shown earlier for CaFe2As2 |133) . obtaining a 
spin- wave fit to the dispersion near the zone bound- 
ary along the [0,1] direction requires very differ- 
ent exchange parameters between nearest-neighbors 
along the spin direction and perpendicular to it. 
Similar dispersions have been measured in the 
isostructural compound SrFe2As2 by Ewings et al. 
jll9j . as shown in Fig. [sj In this figure, results 
are also shown for the paramagnetic state, which 



are quite similar to the ordered state. Two dif- 
ferent approaches to modeling the data are pre- 
sented. Clearly, if one uses a spin-wave model 
that is based on the symmetry of the magnetic 
state, with isotropic nearest-neighbor exchange in 
the paramagnetic phase, then the zone-boundary 
excitations cannot be properly described in the lat- 
ter phase. In contrast, a calculation for itinerant 
electrons within the random-phase approximation 
[117) gets the zone-boundary excitations correctly, 
but has some discrepancy with the temperature de- 
pendence at lower energies. 

Several groups have shown that the anisotropic 
nearest-neighbor exchange can occur due to break- 
ing of the degeneracy between dxz and dyz orbitals 
|134H136] . This orbital ordering [137] is associ- 
ated with the lowering of the lattice symmetry from 
tetragonal on cooling [H [5S] , and it has been ver- 
ified by an angle-resolved photoemission study of 
Ba(Fei_a;Co2;)2As2 |138j . The structural transition 
and its relation to both the antiferromagnetic or- 
der and the superconductivity are of great inter- 
est and have been the subjects of intense study 
[inSl ma [IMISS. two general trends of the 
phase diagrams were established: (i) unless there 
is a first order magneto-structural transition, the 
lattice distortion usually occurs at a higher tem- 
perature {Tg) than the magnetic ordering (Tjv), 
Tg > T/v , and (ii) both Ts and T/v are reduced upon 
chemical substitution, so that both orders tend to 
disappear as the superconducting state develops. 
The orbital ordering presumably occurs at Tg. Ex- 
perimental studies provide evidence for electronic 
anisotropy, often called nematic order, even in the 
tetragonal phase at temperatures not too far above 
Tg [153, 154J. If the nematic order corresponds to 
correlated orbital occupancy, then the anisotropic 
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Figure 8: (a)-(c) Spectrum of magnetic excitations in SrFe2As2 at 6 K, 212 K and 300 K. (d) and (e) show calculations using 
the spin wave theory for local spins with anisotropic, Jia 7^ ■^161 and equal nearest-neighbor coupling, respectively, (f) and (g) 
are the results of the five band itinerant model calculation for the ordered and paramagnetic phase, respectively. Reprinted 
figure with permission from Ewings et al. |119| . © 2011 American Physical Society. 



magnetic exchange might still be applicable. Al- 
ternatively, Wysocki et al. |155) have shown that 
similar effects can be modeled in an effective spin 
model by going beyond Heisenberg couplings (of the 
form JSi ■ Sj) to include a biquadratic term (of the 
form -K{S, ■ S^)^). 

A very useful way of experimentally investigat- 
ing the interplay of magnetic and structural tran- 
sitions in a large single crystal sample used for 
neutron-spectroscopic study on the direct-geometry 
time-of-flight instrument is provided by deploying a 
quasi-Laue technique |112| . In this mode, a quasi- 
white neutron beam with a broad band of incident 
neutron energies centered around some incident en- 
ergy Ei is selected by the pre-monochromating Tq 
chopper. The crystal is aligned with respect to 
the incident beam direction so that Bragg reflec- 
tions of interest appear on the detector. The de- 
tector signal is dominated by the elastic processes 
(diffraction), where the scattering angle is deter- 
mined by the incident neutron energy (or wave- 
length, Xi) and the d-spacing of the set of crys- 
tal planes involved in reflection, in accordance with 
Bragg's law, Xi = 2d sin 9. Such a measurement 
is particularly well suited for studying the relative 
temperature evolution of structural and magnetic 



scattering, which are both present in the diffraction 
pattern at each T. Figure |9] shows an example of 
such measurements performed with Ei w 300 meV 

m- 

The bipartite antiferromagnetic ordering ob- 
served in iron pnictide families is usually described 
as the spin-density wave (SDW) associated with the 
Fermi surface nesting of itinerant electrons [1561 - 
I160j . which corresponds to the {h,k) = (1/2,1/2) 
position in the a&-plane PlI ^nMHTUS] . Such in- 
terpretation is based on the fact that the nesting 
wave vector matches that of the AF order, and 
on the small values of the observed ordered mo- 
ment, which range from « 0.9/is for NdFeAsO and 
BaFezAsa [Mllini], to « 0.4^5 for LaOFeAs [102], 
and to only about « 0.1 /is in Nai-^FeAs |105] . 
Interestingly, the DFT band structure calculations 
typically overestimate the SDW order, predicting 
ordered moment consistent with nearly full spin po- 
larization of an entire electronic band [156] . 

On the other hand, the parent chalcogenide, 
Fei+yTe system, shows a "bi-collinear" magnetic 
ordering with the propagation vector (h, k) — 
(1/2,0) in PA/nmm reciprocal lattice units, which 
does not satisfy the nesting condition. The ground- 
state ordered moment, (/i) < 2/zb |148H152] . al- 
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Figure 9: Maps of the scattered neutron intensity from 
Fei.iTe sample on the ARCS detector bank in the Laue 
mode as a function of the scattering direction at T = 9 K 
[(a) and (c)] and T = 80 K [(b) and (d)]. The top row shows 
structural scattering, while the bottom is magnetic scatter- 
ing. From Zaliznyak et al. |112l . © 2012 American Physical 
Society. 



though larger than in parent fcrropnictides, is sig- 
nificantly smaller than the fully saturated value of 
qUbS (S = 3/2) corresponding to the fluctuating 
paramagnetic moment /Lteff « obtained from 
the Curie- Weiss behavior above 100 K [mi ITT^ . 
In the local-spin picture, the low values of the AF 
ordered moment could arise from frustration and 
low-dimensionality of magnetic exchange interac- 
tions expected in these systems. 

The magnetic form factor provides another im- 
portant clue as to the nature of magnetism in these 
materials. It was recently studied by magnetic neu- 
tron diffraction in the antiferromagnetic SrFe2As2 
[1611 1162j and by the polarized neutron diffraction 
for both the AF BaFe2As2 |163j and the super- 
conducting Ba(Fei_a;Co2;)2As2 |1641 1165] composi- 
tions. The general features revealed by these stud- 
ies are that (i) the magnetization density is only 
weakly anisotropic in wave vector, indicating that 
multiple d-orbitals of the iron atom contribute to 
the observed magnetic moments, and (ii) the Q- 
dependence of the form factor agrees with that for 
Fe atomic orbitals, indicating little hybridization of 
magnetic d— electrons with anions. 

The full value of the fluctuating magnetic mo- 
ment per Fe can be obtained by properly normal- 
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Figure 10: Energy dependence of the imaginary part of the 
dynamical magnetic susceptibility x"(Q,i?), near its maxi- 
mum at {h,k) ^ (0.5,0), in Fei.iTe at 10 K, 80 K and 300 
K, from Zaliznyak et al. 11291 . © 2011 American Physical 
Society. 



izing the spectral intensity of magnetic excitations 
measured in neutron experiment, and applying the 
sum rule for S{Q,E). Zaliznyak et al. |129| car- 
ried out such a program for Fei.iTe and found 
that the low-energy part of magnetic fluctuations 
spectrum, for E < 30 meV, at 300 K already ac- 
counts for the local moment ii^s ~ 3.6/iB. This is 
nearly consistent with fics ~ 3.87/iB corresponding 
to S* = 3/2 with g « 2, which is implicated in the 
Curie- Weiss behavior of the uniform magnetic sus- 
ceptibility. The total neutron intensity, however, 
decreases roughly by a factor 2 upon cooling below 
« 80 K, corresponding to a change in the local spin 
value from S' = 3/2toS'=l. This could be viewed 
as an indication of the Kondo-like screening of the 
local spins by itinerant conduction electrons. The 
spectrum of x"{Qi E) in Fci.iTe near its maximum 
at Q = {h,k) w (0.5,0) at several temperatures is 
shown in Figure |10[ 

^.2. Magnetic correlations in the superconducting 
phase 

Magnetic order in the parent compound gradu- 
ally diminishes with chemical doping, and super- 
conductivity starts to emerge. In some systems 
the two phases are well separated in the phase di- 
agram [140,, .166J . while in others they may coex- 
ist |167H169] . either on an atomic scale or through 
mesoscopic phase segregation. Eventually the AF 
phase is entirely suppressed and the system be- 
comes a good superconductor at low temperature 
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Figure 11: Constant-energy plots of the magnetic excitations in Fei+yTei_i,Sea; projected onto the H-K plane, (a) - (d) are 
data taken from the x = 0.27 non-superconducting sample; (e) - (h) are data taken from the x = 0.49 superconducting sample, 
(i) to (p) are model calculations based on the Sato-Maki function. Reprinted by permission from Macmillan Publishers Ltd: 
Lumsden et al. fT70\ . © 2010 



below Tc. 

For most of their energy band width, the spin 
fluctuations appear not to be strongly affected 
by superconductivity. An example is shown in 
Fig. [TT] where constant energy slices of mag- 
netic scattering intensity from Fei+yTei-^jSex sam- 
ples (a; = 0.27, non-superconducting; and x = 
0.49, superconducting) are plotted. The high en- 
ergy excitations from the superconducting sample 
are qualitatively the same to those observed in the 
non-superconducting sample. A phenomenological 
model based on the Sato-Maki function |171| . a 
form previously applied to studies of critical mag- 
netic scattering in Cr |172j and La2-a;Sr2;Cu04 [77], 
has been used to fit the data here, with parameters 



taken to be energy dependent. 

A similar situation has been observed for high 
energy magnetic excitations in the 122 system. 
Measurements on both the parent compound 
BaFezAsa [HS] (or CaFcaAsa [HTl [SS]) and su- 
perconducting Ba(Fei_^Coa;)2As2 |173[ I174j re- 
veal similar magnetic dispersions. The simi- 
larity between superconducting and parent com- 
pounds is especially clear in the comparison of Q- 
integrated spectral weight between BaFeaAsa and 
BaFe1.9Nio.1As2 as obtained by Liu et al. |125| and 
displayed in Fig. 12 As one can see, the spec- 
tral weight for spin excitations above 100 meV is 
virtually identical. The changes associated with 
doping and superconductivity are restricted to rel- 
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Figure 12; Energy dependence of the imaginary part of the 
dynamical local magnetic susceptibility x" ('^) for BaFe2As2, 
and for BaFei gNio.iAs2 below (filled red circles) and above 
(open red circles) Tc. Reprinted by permission from Macmil- 
lan Publishers Ltd: Liu et al. p^ . © 2012. 



atively smaU energies. The integrated magnetic 
weight, however, does vary substantially among dif- 
ferent Fe-based families. In this case, the sum- 
rule weight corresponds to ^off ~ |125l I126j . 
which is much smaller than in iron telluride. Within 
the experimental error this roughly corresponds to 
S = 1/2, or one unpaired electron per Fe. 

The temperature-dependent changes to the 
magnetic excitation spectra associated with 
the development of superconductivity are sim- 
ilar to those in the high-Tc cuprates. In the 
superconducting phase, a "spin resonance" 
and a spin gap have been observed in vari- 
ous systems, such as the "122" compounds 
Bai_,K,Fe2As2 [173 HZSj, Ba(Fei_,Co,)2As2 
[1771 [T75] . Ba(Fei_:,Ni^)2As2 1^, and the 
"11" compound Fei+^Tei^^Se^r [HH HIS- In 
general the resonance energy Er scales with Tc as 
Er « 4.3 ~ 5.3fcTc [173 [T771 dST]; however this 
scaling appears to break down when pressure is 
applied |183j . The weak-coupling interpretation of 
the resonance and spin-gap in Fe-based systems is 
similar to that applied to the high- Tc cuprates. The 
existence of a resonance has been taken as evidence 
that the sign of the superconducting gap changes 
between Fermi surface pockets that are separated 
by the characteristic antiferromagnetic wave vector 



[1841 I185j . although alternative perspectives have 
been advocated }186j . 

In both the 122 and 11 systems, the resonance 
always occurs around the in-plane wave vector 
(0.5,0.5) despite different ordered ground states in 
their parent compounds. The excitations at higher 
energies appear to disperse only in the direction 
transverse to this wave vector. It has been pro- 



the presence of orbital correlations [187j ; note that 
it has also been proposed |l88j that orbital fluctu- 
ations, rather than spin fluctuations, might medi- 
ate superconductivity. In contrast. Castellan et al. 
|176| have observed a longitudinal splitting of the 
spin resonance in Bai_2,K^Fe2As2 at high doping. 
Theoretical RPA-type calculations taking account 
of imperfect Fermi surface nesting can reproduce 
this effect [T75|. 

Both the intensity and energy of the resonance 
mode are found to be affected by temperature |178j 
and external magnetic field |189|. I190j , demonstrat- 
ing an intimate correlation between the supercon- 
ducting electron pairing and magnetic excitations 
in the Fe-based superconductors. Nevertheless, the 
microscopic origin of the resonance mode is still not 
fully understood. Is it a singlet to triplet excita- 
tion arising from quasi-particle scattering, or sim- 
ply a modification of existing magnetic excitations 
by the establishment of superconductivity? If the 
former is the Zeeman splitting of the triplet 

mode is expected with the application of an exter- 
nal magnetic field. In the Fe-based superconduc- 
tors, because of the reduced compared to the 
cuprates, the resonance energy E^ is also consid- 
erably lower, making it easier to pursue this prob- 
lem. A number of inelastic neutron scattering stud- 
ies have been performed, aiming at understanding 
the change of the resonance mode under a magnetic 
field [rSnfTSMTOT] . but so far there is no conclu- 
sive evidence on a possible Zeeman splitting and 
the question of the origin of the resonance remains 
open. 

As mentioned previously, orbital order involv- 
ing broken degeneracy of the d^z and dy^ states 
has been predicted |134H136] and experimentally 
confirmed [138] in the lowered-symmetry structural 
phase of iron pnictides and chalcogenides. Suffi- 
cient chemical substitution leads to a restoration of 
tetragonal lattice symmetry, and hence orbital or- 
der should be suppressed; nevertheless, freezing of 
local orbital correlations may still occur |135l 1193] . 
Although such effects are difficult to detect di- 
rectly, they should have an impact on the mag- 
netic response. Indeed a dramatic change in the 
dispersion of low energy magnetic excitations has 
been observed in the superconducting 11 compound 
[1921 1194] ■ at a temperature of around 3Tc (see 
Fig. 13). At about the same temperature scale, 
an abnormal in-plane expansion [152, 192J is ob- 
served, which would be consistent with freezing of 
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Figure 13: (a) to (c) Low energy magnetic excitations in a 
Feo.96Nio.04Teo.5Seo. 5 superconducting (Tq = 8 K) sample, 
(d) is the same data measured from a non-superconducting 
Feo.9Cuo.1Teo.5Seo. 5 sample, taken at 4 K; from Xu et 
al. [T92] . 

local orbital correlations. In addition, thermoelec- 
tric power [195] , and optical conductivity |1961 1197) 
measurements both indicate a change of electronic 
states near the Fermi level in this temperature 
range. Although the exact role of orbital corre- 
lations in the magnetic and superconducting corre- 
lations are not fully understood, this is a direction 
that deserves further attention. 

Finally, we remind the reader that this is far from 
an exhaustive review of the field. In particular, we 
note that there are some newer systems, such as 
(Li,Na)FeAs and A^Fea-^Ses (A = K, Rb, Cs, Tl), 
that we have neglected for the sake of brevity. For 
pointers to some of the neutron work on these very 
interesting materials, the reader may wish to con- 
sult the reviews [32] and |5T] . 
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